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We consider the probability by which quantum phase measurements of a given precision can be 
done successfully. The least upper bound of this probability is derived and the associated optimal 
state vectors are determined. The probability bound represents an unique and continuous transition 
between macroscopic and microscopic measurement precisions. 
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The classical picture for the evolution of a single-mode 
electromagnetic field is simply determined by an ampli- 
tude (specifying the strength of the field) and a phase 
(specifying the zeros of the field). On the other hand, the 
concept of electromagnetic phase as an observable quan- 
tity is a long-standing problem of quantum optics and 
it has been the question whether there exists a phase 
observable that is canonically conjugate to the number 
observable for a single-mode field. The quantum mechan- 
ical description of phase was first considered by London 
[H and Dirac . An obvious way of defining an operator 
for the phase is by polar decomposition of the photon 

annihilation operator a = e^vfi. The phase operator 
4> defined in this way is equivalent to that considered 
by Dirac jjj, who obtained the commutator [cj>, N] — i by 
employing the correspondence between commutators and 
classical Poisson bracket. Formally, this would imply the 
uncertainty relation 
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with (70 and ct/v are the standard deviations of </> and N. 
The difficulties of Dirac's approach were clearly pointed 
out by Susskind and Glogower [3J. Firstly, the relation 
([!]) would imply that a well-defined number state would 
have a phase standard deviation greater than 2ir. This is 
a consequence of the fact that Dirac's commutator does 
not take account of the periodic nature of the phase. Fur- 
thermore, the exponential operator e 1 ^ derived from this 
approach is not unitary and thus does not define a Her- 
mitian operator. This is why it is often accepted that 
a well-behaved Hermitian phase operator does not exist 
0, 0]- Therefore, arguments based on the Heisenberg 
relation {J) cannot hold in general. 

Actually, the standard deviation offers a reasonable 
measure of the spread of values when the distribution 
in question is of a simple "single hump" type. In par- 
ticular it is a very good characteristic for a Gaussian 
distribution since it measures directly the half-width of 
this distribution. However, when the distribution is not 
of a simple type (for example, has more than one hump) 
the standard deviation loses much of its usefulness as a 
measure of uncertainty. 

The aim of the present contribution is to introduce the 
probability by which a successful phase measurements of 



a given precision can be done. The least upper bound 
of this probability is determined and the corresponding 
(optimal) state vectors are computed. 

In order to specify phase measurements, the probabil- 
ity distribution for the measurement result can be de- 
termined using positive operator-valued measures. This 
approach was first considered by Helstrom |j§], and is 
also considered in @, 0|- Precisely, let TL be a com- 
plex separable Hilbert space, (|n)) n >o an orthonormal 
basis, and N = J2n°=o n \ n )( n \ the associated number 
observable. If the phase density F$ treats all phases 
equally, it should be invariant under phase translation, 
e N F^e — F^+g, generated by the number observ- 
able. In this case, the general form of F^ is 
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where (c nm ) is the associated phase matrix. For the in- 
tegral of the probability to equal 1, we must have 



F rh dcb = l. 
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Applying this to ((2j) above we find that 

oo 

^c„„ \n)(n\ = 1. (4) 
n=0 

This means that the diagonal elements c nn must all be 
equal to 1. The additional condition of positive definite 
probabilities, together with the above result means that 
all of the c nm must have absolute values between and 
1. In general, real measurements will give smaller val- 
ues of c nm , and the closer these are to 1 the better the 
phase measurement is. In 0| it is shown that the addi- 
tional condition that a number shifter does not alter the 
phase distribution gives c nm — 1, corresponding to the 
canonical measure [9|| 
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An alternative derivation of this result is by using the 
maximum likelihood approach jf|. Note that <[5j) may be 
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expressed by dE^ = \4>){4>\ d<fi, where 
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With reference to 10, we now define the precision 
Aa £ [0, 2tt) of a phase measurement corresponding 
to the vicinity A a — [a — 4p, a + As) of any value 
a € [— 7r,7r). The probability of a phase measurement 
with 4> £ A a , made on a state described by a density 
operator p, is given by 



P a {Aa) = Tr[pE,p{A a ) 
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On the other hand, the probability to measure a pho- 
ton number n g B k is given by 



where 



P fc (Afc) = Tt[ P E N (B k )} 
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is the value of the spectral measure En on a set B k of 
positive integers. In order to introduce the precision 
Afc e N by which the photon number is measured, we 
note that the photon number is bounded from below. 
Therefore, we define the right-sided vicinity B k C N of 
fc by B k — {fc, k + 1, k + Afc}. In this definition, the 
integer fc is the smallest element of Bk . Alternative def- 
initions are also possible but typically lead to necessary 
readjustments in certain cases for fc < Afc. However, 
it will be seen later that our results are not dependent 
on the particular subscription. For technical purposes we 
apply the definition of the minimum integer-subscription. 
In the case of pure states p = \ip) (tp\, we obtain the prob- 
ability 
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and ip n = (n\ip) is the number-space amplitude of ip. 

Now, we consider the case with an initial photon num- 
ber preparation of a state ip, A single mode is supposed 
to emerge in a state according to 
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Afterwards, the number of photons is given with precision 
Afc. In this situation the uncertainty principle suggests 
that the more accurately the number is measured the 
greater is the perturbation of the phase of the outgoing 
state. The conditional probability V a k(Aa \ Afc; ip) to 
measure phase <p € A a , on the state transformed by the 
initial number measurement, is given by 



V ak (Aa\ Afc;V) 
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We now ask for the least upper bound Xq of the mea- 
surement probability lfT2|) and we end up in a variation 
problem in Hilbert space with three degrees of freedom. 
For fixed precisions Afc and Aa we are searching for the 
supremum of I112II by variation of the parameters a, k and 
the state vector ip of the photon. Actually this variation 
problem is translation and rotation invariant in Hilbert 
space and we can simply chose fc = and a = with- 
out loss of generality. After all, we have to consider the 
following expression 



A = sup V a(Aa\Ak;ip) 

V>6H\{0} 
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and by using J5]) and J9]) we explicitly obtain to following 
expression 
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Applying the Cauchy-Bunyakovsky inequality we obtain 
the following general upper bound 
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for every Aa G [— n, ir) and integer Afc > 0. In fig. 1, the 
set of impossible measurement processes is expressed by 
the grey shaded triangle. In order to reach even tighter 
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Figure 1: Maximum eigenvalue Ao for Afc = 0, 1, 2, 3, oo (from 
left to right). Phase measuring processes with probabilities 
larger than Ao do not exist. 

bounds we explicitly computed the integral in (fl4|) and 
applied certain trigonometric identities to obtain the fol- 
lowing expression 



Poo(Aa | Afc; tp) = ^ G nm ipni>* m 
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with normalization condition X)n=o \^n\ 2 = 1 and kernel 
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for < n,m < Afc, otherwise. Obviously, G is self- 
adjoint and positive definite. According to lfl"3|) and (fl6|l . 
the least upper bound is given by the operator norm of 
G, i.e. 



Ao = \\G\\, 
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and this norm is identical to the largest eigenvalue of G. 
In order to obtain the eigenvalues of G we have to solve 
the following linear equation for n = 0, 1, Afc 
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for s = 0,1,..., Afc. This type of eigenvalue problem 
has been extensively discussed in |ll| (see also references 
therein). All eigenvalues are distinct, positive and may 
be ordered so that 1 > Ao > Ai > ... > Aa/c- In the 
non-trivial case of Afc > we computed A numerically. 
For increasing values of Afc, the corresponding bounds 
approach very fast to the asymptotic case Afc — > oo, see 
fig. 1 (right most continuous line). For the computation 
of the asymptotic case we introduced the equidistant de- 



, Tn = 0, 1, Afc with increment 
After substitution into lfl9l) and a few 



composition q m ■ 
Sq m = %n+i — Qn 
algebraic manipulations, the discrete eigenvalue problem 
approaches to the following homogeneous Fredholm inte- 
gral equation of the first kind 

- 1 sin(§ f(z - z')) 
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with \z\ < 1 and 
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From standard theory we know that there are solutions 
in L 2 ({— 1, 1]) only for a discrete set of eigenvalues, say 
Ao > Ai > ... and that as v — > oo, \ v — > 0. It 
should be noted that the eigenvalues explicitly depend 
on the parameter £ and corresponding to each eigen- 
value there is a unique (up to normalization) solution 
<p( v '{z) — S'o ! y(7r£/2 J _z) called angular prolate spheroidal 
wave function UJ, III! ■ They are continuous functions of 
£, for £ > 0, and are orthogonal in (—1, 1). Moreover, they 
are complete in L 2 {[— 1, 1]). The corresponding eigenval- 
ues are related to a second set of functions called radial 
prolate spheroidal functions, which differ from the angu- 
lar functions only by a real scale factor [ll|. Applying 



the notation of [lO, [ll| , these eigenvalues are 
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with v — 0, 1,2, ... The properties of the discrete eigen- 
value spectrum is discussed in [l3|. Here, we are mainly 
interested in the properties of the largest eigenvalue 
A (£)- It is monotonically increasing and approaches 1 
exponentially in £. For small values of £ there is the 
asymptotic behavior Aq(£) ~ £. 
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